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Abstract 


Introduction 

As light passes through a optical system the reflections and refractions will in general change the polarization state 
of the light. If we assume that all of the materials in the thin film coatings and substrate are isotropic and 
homogeneous then calculating the amount of "instrumental ' polarization is a relatively straight forward task. In the 
following sections we will present all of the steps required to perform a "polarization ray trace" calculation for a single 
ray and monochromatic and hence polarized light. The thin film portion of the calculation is also shown. The reason 
for explicitly showing the thin film equations is that there are sign conventions imposed on the boundary value 
equations by the orientation and handedness of the various coordinate frames which are attached to the geometric rays. 

The attenuation of light through a optical system, is relatively simple, and requires at the very least a lens (average) 
reflectivity or transmissivity. Determining the polarization sensitivity of a optical system is still relatively straight 
forward requiring at least a knowledge of the behavior of the “s” and “p” components at each interface for the chief 
ray. Determining the thin film induced aberrations of a optical system are somewhat more demanding. Questions 
about the arithmetic sign of the phase factors and how this relates to the overall “OPD” of a ray are ubiquitous. Many 
rays are required to construct a wavefront. Thin film codes which modify the OPD s of rays are a requirement for 
this last mentioned computation. This requires a consistent scheme of coordinate frames and sign conventions and 
is probably the most demanding task of a polarization ray trace. 

Only the electric field will used in the discussion. This is not a restriction as the Stokes parameters are functions of 
the electric field. The following does not attempt to explain, but only to present all of the required concepts and 
formulas. 


Maxwell's Equations in a Conducting Medium 
It is very convient simply to state that in an isotropic medium with 


we have 


e = dielectric constant 
fx = permeability and 
a = conductivity 
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Introduction 

As light passes through a optical system the reflections and 
refractions will in general change the polarization state of the 
light. If we assume that all of the materials in the thin film 
coating and substrate are isotropic and homogeneous then 
calculating the amount of "instrumental ,, polarization is a 
relatively straight forward task. In the following sections we 
wi.li present all of the steps required to perform a "polarization 
ray trace" calculation including the thin film portion of the 
calculation. The reason for explicitly showing the thin film 
equations is that there are sign conventions imposed on the 
boundary value equations by the orientation and handedness of the 
various coordinate frames which are attached to the geometric rays. 

The following does not attempt to explain, but only to present all 
of the required concepts and formulas. h , *-■ , /_ 

Maxwell's Equations in a Conducting Madium 

In an isotropic medium with 

e = dielectric constant 
H = permeability and 
a = conductivity 


then 
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Wave Equation 

If the electric field is given by 
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and the row vector by 


g- = (e; e;) = (B| 

with E* denoting complex conjugation and E + Hermitian conjugation. 

Local Rav Coordinate Frames 



Figure 2: Local Coordinate Frames. 


At each interface we have a incident ray , a reflected ray ic r and 
a transmitted ray and a surface normal fl , which is oriented such 

that the scalar product > 0 . The local s and p directions for 
each of these rays is given by 
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where is perpendicular to the plane of incidence and tangent to 
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or in a more compact notation simply by 
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where the subscript R 2 indicates the incident surface. 

Boundary Conditions 



Figure 4: Two interfaces. 


At each interface, within a film stack, the general form of the 
equations expressing continuity of the tangential components of the 
E and H fields is similar to the equations at the first interface. 

{E x + E'^ x ft = ( E 2 + e! 2 ) x ft ( 16 ) 
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Comparing equations (19 & 20) and (24 & 25) we see that at each 
interface we have equations of the form 
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With the following substitutions for the s-component 

Z= E 
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C = — COS0 
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and for the p-component 
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Multi-Laver Stack Equations 

Referring to Figure (5) we have the following equations (for either 
the s or p components) at the first interface 
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At Succeeding (1,1+1) interfaces we have, with 
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we have 



Matrix form of the Equations 

Equations (32-38) can be rewritten as 



or with A and B representing 2 by 2 matrices and "ket" 
appropriate column vector we have equations (42) 
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The form of the inverse matrices is 
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Matrix Form 


We can express the reflected and transmitted fields after the m th 
surface as 
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And similarly to equ.(15) and (15.1) this can be written as 
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Polarization Ray Trace 

From equations (15.1) and (55) we see that at each surface the 
expression 

= S^pn-l) (56) 

connects the electric field before and after each surface. The 
expression 

pn+l) ~ (57) 

shows how the incoming plane wave electric field is modified by the 
intervening surfaces. 


Energy Flux 


The magnitude of the energy flux, J m in a pencil of light, after 
the m* surface is given by (to within a constant of 
proportionality) 
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and from this we see that the maximum and minimum occur at 

0 and 6 + it/ 2. 


Phase Effects 


The transmitted intensity is given by 
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where s and p are complex field quantities. The Z matrix is 
Hermit ian and is the product of the two matrices in equ. (60) where 
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from equ. (57) . Now 
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the t denoting transpose and the * denoting complex conjugation. 
The general forms of the matrices are 
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hence 
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and the time average is given by 
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The real part of the time integral is 
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From this we see that only when 


tCw-ojO < -5. 


is the integral different from zero. For 


1=1 micron 


a band pass filter of AX - 10 nanometer we get 
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